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Abstract 


The  problem  of  estimating  the  state  of  a  target  in  the  presence  of 
measurements  of  uncertain  origin  has  received  a  great  deal  of  attention 
recently.  When  this  origin  question  cannot  be  resolved  with  certainty, 
one  can  only  make  probabilistic  inferences  as  to  which  detection  (and 
thus  measurement)  originated  from  the  target  of  interest.  The  investi¬ 
gation  reported  here  deals  with  a  new  class  of  problems  characterized 
by  the  following:  the  correct  measurement  arrival  (detection)  times 
for  a  target  of  interest  occur  according  to  a  stochastic  process. 
Another  stochastic  process  governs  the  arrival  times  of  the  false 
alarms.  Thus  detections  occur  one  at  a  <lme  and  while  some  of  them 
can  be  discarded  as  not  having  originated  from  the  target,  the  remain¬ 
ing  ones  cannot  be  associated  with  certainty  with  the  target.  This 
problem  is  motivated  by  the  fact  that  in  some  tracking  problems 
detections  from  the  target  of  interest  occurs  on  an  irregular  basis. 

A  procedure  that  associates  probabilistically  these  measurements  to 
the  target  is  developed  together  with  a  corresponding  estimator. 

The  optimal  estimator  as  well  as  a  number  of  suboptlmal  algorithms  that 
are  real  time  lmplementable  are  presented  together  with  simulation 
results.  The  simulations  also  indicate  that  the  probabilistic  data 
association  filter  with  t lme-of-arrival  information  is  significantly 
superior  to  the  filter  which  uses  only  measurement  location  information 
for  probabilistic  data  association. 

The  second  part  deals  with  the  problem  of  Joint  probabilistic 
data  association.  In  this  case  there  are  several  targets  and  each 
measurement  could  have  originated  from  one  of  these  targets  or  none. 

The  approach  to  this  problem  is  to  obtain  joint  probabilities  for  all 
the  feasible  associations  of  measurements  to  targets.  An  example  is 
provided  that  illustrates  that  the  results  are  substantially  different 
from  the  case  when  marginal  probabilities  are  obtained  directly. 
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0 .  Introduction 

This  report  is  organized  In  two  parts.  The  first  part  deals  with  the 
case  of  a  single  target  from  which  measurements  are  obtained  at  times  governed 
by  a  stochastic  process,  (i.e.  asynchronous  measurements)  with  false  measure¬ 
ments  arriving  at  times  governed  by  another  stochastic  process.  The  purpose 
of  this  study  was  to  develop  the  methodology  of  incorporating  the  time-of- 
arrival  Information  into  the  data  association  procedure.  The  results  indicate 
both  the  feasibility  as  well  as  the  significant  performance  Improvements  of 
the  probabilistic  data  association  method  with  t ime-of-arr lval  information. 

The  second  part  deals  with  the  probabilistic  association  of  an  arbitrary 
number  of  measurements  to  an  arbitrary  number  of  targets.  The  result,  called 
Joint  Probabilistic  Data  Association,  accounts  for  all  the  possible  origins  of 
each  measurement.  References  are  presented  following  each  part. 

This  study  has  been  motivated  by  several  practical  problems  encountered 
mainly  In  passive  sonar  tracking.  The  algorithms  developed  In  the  course  of 
the  study  are  feasible  for  real-time  Implementation  and  will  allow  more  efficient 
use  of  the  available  Information  In  tracking  systems. 
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Part  I 


TRACKING  WITH  MEASUREMENTS  OF  UNCERTAIN  ORIGIN 
AND  RANDOM  ARRIVAL  TIMES 


T 


1 .  Int  roduet Ion 


The  problem  of  estimating  the  state  of  a  target  In  the  presence  of  measure¬ 
ments  of  uncertain  origin  has  received  a  great  deal  of  attention  recently.  This 
attention  was  motivated  by  the  fact  that  when  tracking  with  certain  sensors  one 
can  observe  several  detections  whose  origin  might  not  have  been  the  target  of  interest. 
When  this  origin  question  cannot  be  resolved  with  certainty,  one  can  only  make 
probabilistic  inferences  as  to  which  detection  (and  thus  measurement)  originated 
from  the  target  of  interest.  This  problem  arises  mainly  In  passive  sensors  of 
acoustical  or  optical  type  and  can  be  due  to  clutter  or  false  alarms  or  inter¬ 
ference  from  undesirable  targets.  In  an  environment  where  there  are  several 
targets  as  well  as  false  alarms  the  ambiguities  in  the  association  of  the 
measurements  to  the  targets  can  become  a  rather  severe  problem. 

The  previous  work  related  to  this  multi-target  tracking  problem,  where  there 
ran  be  targets  of  interest,  some  of  no  interest  and  also  false  alarms  (false 
targets),  was  recently  surveyed  in  |B1].  The  approaches  discussed  in  the 
literature  were  either  of  the  likelihood  function  or  Bayesian  type.  They 
assumed  that  at  predetermined  times  "snapshots"  of  the  environment  were  obtained 
with  several  detections  (and  the  corresponding  measurements)  that  were  to  he 
associated  into  tracks. 

The  investigation  reported  here  deals  with  a  new  class  of  problems  charac¬ 
terized  bv  the  following:  the  correct  measurement  arrival  (detection)  times  for 
a  target  of  interest  occur  according  to  a  stochastic  process.  Another  stochastic 
process  governs  the  arrival  times  of  the  false  alarms.  Thus  detections  occur  one 
at  a  time  and  while  some  of  them  can  be  discarded  as  not  having  originated  from  the 
target,  the  remaining  ones  cannot  be  associated  with  certainty  with  the  target. 

A  procedure  that  associates  probabilistically  these  measurements  to  the  target 
is  developed  together  with  a  corresponding  estimator.  This  problem  is  motivated 
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by  the  fact  that  in  some  tracking  problems  contact  with  the  target  of  Interest 
occurs  on  an  irregular  basis.  A  related  problem  was  studied  by  Nahi  |N1]  and 
.laffer  and  Gupta  (J1,J2)  where  they  considered  the  measurements,  obtained  at 
fixed  times,  as  being  either  a  linear  combination  of  the  states  plus  noise  or 
(zero-mean)  noise  alone.  An  estimation  problem  with  measurements  occuring  at 
random  times  was  considered  in  Snyder  and  Fishman  ( S2  ]  .  A  combined  estimation  and 
control  problem  with  some  of  the  measurements  occuring  at  random  times  was 
studied  in  Rhodes  and  Snyder  (Rl).  However,  in  these  investigations  there  was 
no  Origin  uncertainty  about  the  observations. 

The  problem,  formulated  in  Section  2,  models  the  arrival  times  of  the 
measurements  as  arbltrarv  renewal  processes.  The  full  memory  Probability  Data 
Association  procedure,  which  carries  all  the  origin  hypotheses  up  to  the 
current  time,  is  developed  for  this  model  In  Section  3,  and  the  corresponding 
estimator  in  Section  4.  Suboptlmal  procedures  are  presented  in  Section  !>. 

The  simplest  of  them,  the  Combined  Single  Hypothesis  (CSH)  filter,  is  only 
slightly  more  expensive  to  Implement  than  the  Standard  Kalman  Filter  (SKF). 

The  simulation  results  presented  in  Section  A  compare  the  SKF,  CSH  and  a 
limited  memory  version  of  the  PDA  filter  that  retains  the  N  most  likely 
hypotheses.  A  significant  Improvement  in  performance  has  been  observed  for 
the  CSH  over  the  SKF.  The  additional  improvement  obtained  over  the  CSH 
when  carrying  N“8  hypotheses  was  not  so  dramatic. 
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2.  Formulation  of  the  Problem 

It  Is  assumed  that  there  is  one  target  of  interest,  whose  state  x,  an  n- 
vector,  evolves  in  continuous  time  according  to  the  equation 


x ( t )  -  A  x(t)  +  B  w(t) 


(2.1) 


where  A  and  B  are  known  matrices  and  w(t)  is  a  zero-mean,  white,  Gaussian  pro¬ 
cess  noise  with  known  covariance 


w ( t )  w'(t)  •  q  6(t- t) 


(2.2) 


A  measurement  corresponding  to  a  detection  that  originated  from  the  target 
being  tracked  is  given  by  the  equation 


z(tfc)  -  H  *(tk)  +  v(t  ) 


(2.3) 


where  t,  is  the  time  at  which  the  detection  occurs  and  v  is  a  zero-mean  white 

k 

Gaussian  noise  sequence  with  covariance 


E  v(.k>  V(tj)  -  «kJ 


(2.4) 


and  independent  of  the  process  noise. 

The  set  of  measurements  and  the  corresponding  arrival  times  are  denoted  as 


z  - 


(2.5) 


T  -  Vi.i 


(2.6) 


The  minimum  mean  square  error  estimate  of  the  state  is  then 


x  ( tfc ! tk>  "  E[x(tk) !zk,  Tk] 


(2.7) 


The  arrival  times  of  the  measurements,  both  correct  and  false  ones  are 
modeled  as  renewal  processes. 


The  difference  between  "sensor  time"  and  "target  time"  has  been  pointed  out  in 
[FI].  Here  we  do  not  make  this  distinction. 
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The  pdf  of  the  arrival  time  of  the  correct  measurements  is 

p(t)-  1  >  '(  1  (2,8) 

where  t  .  f ••  the  time  of  the  previous  correct  detection. 

*•»  i 

The  pdf  of  the  arrival  time  of  the  false  measurements  is 


‘  V'-'i.o’  '  '  Vo 

where  t.  ^  Is  the  time  of  the  last  incorrect  measurement  prior  to  t.  The 
symbols  K  and  will  be  used  for  the  distributions  corresponding  to  the 
densities  in  (2.8)  and  (2.9),  respectively. 


When  a  new  measurement  is  received,  sav,  at  t,  ,  a  validation  region  is 

k 

set  up  around  the  predicted  value  of  the  measurement  at  t^  based  upon  past  data. 

If  the  new  measurement  is  outside  this  region  it  is  discarded;  otherwise  it  will 

he  associated  probab  1 1  1st ica 1 1 v  to  the  track  and  used  to  update  the  target's 

estimate.  This  is  discussed  in  detail  later. 

If  the  measurement  is  false,  e.g.,  from  clutter  or  false  alarm,  the  spatial 

distribution  of  this  measurement  is  assumed  independent  of  past  occurences 

and  uniform  within  the  surve  1 1  lance  region  (with  hvpervolumo  V  ) 

s 

P[Z(tk)!Vo'  zk’''  Tk~']  •  pl2(tk)!\,o]  ■  1/Vs  (2-10) 


where  we  define  as  follows  the  events 


(i  4.  (detection  at  t^  did  not  originate  from  the  target) 
j  ^  {detection  at  t^  originated  from  the  target) 


(2.11) 


The  information  to  be  used  in  calculating  the  probability  of  each  measure- 


ment  being  correct  (the  probabilistic  data  association  -  PDA)  is  their  arrival 
time  and  location.  In  order  to  use  all  the  information,  past  history  hypotheses 

k.l  m  {  xk-l.qf  ^  }  (2.12) 

have  to  be  stored  and  their  probabilities  evaluated.  Each  hypothesis  is  an 
assumption  as  to  the  origins  of  all  the  past  measurements.  This  "full-memory" 

PDA  is  developed  in  the  next  section.  It  is  somewhat  similar  to  the  "full  scan 

* 

back"  approach  of  Singer  et  al. [SI).  The  main  difference  is  that  here  one  deals 
with  random  arrival  times.  Since  the  requirements  of  the  full-memory  PDA  grow 
with  time,  " l imi ted-memorv"  versions  of  it  are  also  presented  later. 


This  terminology  Is  radar  oriented  where  several  detections  are  obtained  at  each 
scan,  l.e.,  at  fixed  times. 
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3.  The  Probabilistic  Data  Association  for  Measurements  with  Random 
Arrival  Times  -  Full  Memory _ _ _ _ _ 

The  conditional  mean  of  the  state  (2.7)  can  be  written  using  the  total 

probability  theorem  as 


where 


X(tk,tk) 


•Ee[«(t k)|xk*1.  zk,  Tk]  s k ’ 1 
i 


6k>1  2  p{ xk ’ 1 ! zk . 


Tk) 


(3.1) 


(3.2) 


is  the  conditional  probability  that  the  i-th  past  history  Is  correct.  The 
evaluation  of  the  probabilities  (3.2)  constitutes  the  PDA  method.  This  method 
associates  past  sequences  of  measurements  to  the  target  In  a  probabilistic 
manner,  hence  Its  name  [B2.S1 ] . 

Each  estimate  of  the  state  vector  conditioned  upon  a  particular  hypothesis 
on  the  r.h.s.  of  (3.1)  Is  obtained  by  a  standard  linear  (Kalman)  filter. 

The  probabilities  (3.2)  will  be  evaluated  next  using  Bayes'  rule  and  the 
model  of  the  measurements  presented  In  the  previous  section. 


k,l  k  k  k- 1 , q  _k-  1  .  _k- 1  , 

?(t  Z  .  T  )  -  P( x  \  Xk  j i Z  .  *(tk),  T  .  tk) 

- 1  ,  ,  ,  i  k- 1 , q  k-1  k-1,  _ ,  k-1 ,  q,k- 1  k-1 . 

"  C  P  ^k, J ’ *  £k  ’  Ck  X  *Z  *T  }  P{x  Z  -T  } 


-  c  1  p[z(tk)ixk  1,q.Xk  j*  Zk  1 ,Tk  1 , tk ) 

•  p,Xk<J.tklxk"l,q.  Zk_1,  Tk_1)  P{xk’1,q|Zk‘1,  Tk_1) 


(3.3) 


where  c  Is  a  normalizing  constant.  Each  expression  on  the  r.h.s.  of  (3.3)  will 
be  evaluated  separately. 


The  first  term  Is  the  pdf  of  the  location  of  the  last  measurement  and  is 

* 


either  normal  (for  a  correct  measurement,  i.e.,  J"l)  or  uniform  (for  an  Incorrect 
detection,  l.e.,  J«0) 

This  follows  from  the  1 lnear-Gauss lan  model  (2.1)— (2.4) . 
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pU(tk)|xk"1,q.Xktj.  Tk_1.  t, 


<*[z(tk);  *Uk  tk  l,  xk"1,q).  S(tk|tk_1.  ^  i 

(  v'1  if  j  -  o 

'  u 


H.q, 


(3.4) 


if  J  -  0 


where  ^(z;  r,  S)  is  the  normal  density  with  variable  z,  mean  z  and  covariance 
matrix  S  (which  is  the  innovation  covariance  in  this  case  conditioned  on  a 
particular  hypothesis,  to  be  discussed  In  more  detail  later);  \  is  the  hyper- 

volume  of  the  surveil lance  region. 

* 

The  second  term  is  the  Joint  probability  of  the  event  *  (that  the  detection 

*  •  J 

.it  t  is  of  tvpe  1)  and  the  arrival  time  t,  oi  this  detection  condit ioned  on  a 
k  k 

past  hlstorv  hypothesis.  Note  that  t  is  the  realization  of  the  random  variable 

which  is  "the  arrival  tira«  of  t  tie  first  detection  •  f  tvpe  J  after  t 
k-l.J  k-1 

Therefore  this  Joint  probability  can  he  written  as 

k- 1  q  k-1  k-1 

»MxktJ.fk  X  *q.z  .T  ) 

k-l,q  k-1  k-1 , 

’  ‘k.J  k-l.J  *  V  *  '  1  ' 


k-l.J 


k-lfq  k-1  k-1 

t ,  *  •  /.  ,  I 

k 


(  i.S) 


where  the  density  above  is  subscripted  to  Indicate  the  random  variable  it  describes. 
Denote  bv  ••  the  "last  arrival  time  of  a  detection  of  tvpe  J  prior  to  t  "  (in 

k.J  K 

k  —  1  *  q 

the  entire  surveillance  region;  this  is  yielded  hv  the  conditioning  event  x  >* 


Then  one  has 


>  it  7.k_1  Tk"' 

T  1  L  ’  ‘  * 

k-l.J 


CJ  fJ(tk  “k.J* 


(J.M 


This  is  a  mixed  probability  of  an  event  and  a  continuous  random  variable.  The 
notation  used  is  the  one  for  events,  i.e.,  J’. 


where  the  notat ions  of  (2.8),  (2.9)  are  used.  The  normalization  constant  c 
Is  needed  because  one  has  a  truncated  density:  t  must  be  larger  than  j. 


A 


J  k  k.j  k  I  k-1  k.J 


(3.7) 


Since  event  x  Is  equivalent  to  the  event  that  the  first  detection  from 

k  •  J 

* 

process  1  after  t  comes  later  than  the  one  from  process  ),  the  lirst  proba¬ 
te-  1 

bllltv  on  the  r.h.s.  of  (J.S)  tan  be  written  as 

k-l,q  k-1  k-1 . 

P,V.J  'k-1 ■  'C  • z  • T  ' 

.  p|,  >  ,  ■,  -  I  ,  7.k‘‘,  Tk''l 

k-1 . J*  k-l.j  k-l.j  k* 


*-Y  “k-fk.j'1 


(3.8) 


where  the  denominator  Is  again  a  normalization  constant  c(*  similar  to  (3.7). 
Note  that,  for  any  x^  one  of  these  two  normalization  constants  will  be  unity 

The  last  term  In  (3.3)  is  available  from  the  previous  update.  This  com¬ 
pletes  the  calculation  of  the  a  posteriori  probabilities  for  each  past  history 
hypothesis . 
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4.  The  hstimator 

Once  the  hypothesis  probabilities  have  been  obtained,  the  m.m.s.e. 
estimate  oi  the  target's  state  results  iron  (1.1).  Kach  hypothesis-conditioned 
est Innate  in  (i.l)  Is  obtained  trom  a  standard  Kalman  filter  (which  is  the 
optimum  estimator  when  conditioning  on  a  hypothesis)  together  with  an  associated 
covariance  matrix.  Note  that  the  conditional  density  of  the  target's  state  Is 
a  Gaussian  mixture  with  ti  '  the  weightings.  The  covariance  corresponding  to 
the  estimate  (3.1)  is  (Sl.Bl) 

>-«kiy  -E*k-‘  px„iv 

»dk  >k.  >k'‘>  t'<tk|tk.  »k'S 


x(t 


k  1 


k>  X’(tk 


V 


Ek.l  L  .  k,  1 
|P(tk  k'  ‘  ) 


*  (x-(tk  V  ‘k,i)  ■  x(tk  V1  lx(tk’V  xK,1)  -  x(tk  V1  i 

(4.1) 


k,  1 , 


.1 


In  order  to  reduce  the  computational  requirements  one  can  igti  <re  in  the 

PDA  measurements  that  are  "far"  trom  tin  predicted  measurement  -  i.e.,  one  can 

use  .i  validation  region  that  eliminates  the  measurements  which  would  be  deemed 

by  the  PDA  as  having  a  very  low  probability  of  being  correct.  Since  to  set  up 

an  exact  validation  region  based  on  the  Gaussian  mixture  or  a  different  one  for 

each  hypothesis  is  not  practical,  an  elliptic  region  can  be  calculated  assuming 

the  state  is  approx lma t e 1 v  normally  distributed  with  covariance  as  in  (4.1).  The 

innovation  covariance  at  t,  is  then 

k 
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S(,k  Vl>  -  H  P«kllk.l>  H  ♦  R„ 


(4.2) 


fhe  state  prediction  covariance  Is  obtained  from 


p<'k'Vl’  ‘  “Wl*  '■“v-rVl'  ‘  “k'Vl'  *  "‘Wi’ 


(4.3) 


where  *  Is  the  transition  matrix  corresponding  to  A  In  (2.1);  P ( t  _ ^ j t k_ ^ )  is 
given  by  the  equivalent  of  (4.1)  at  t  ;  the  covariance  Q  Is  the  contribution 
of  the  continuous  time  process  noise  to  the  state  prediction  covariance 


«V‘k-i’  ‘ 


}(t,  ,t)  R  i)  B  4  ( t  ,  t )  d  t 

k  k 


(4.4) 


Th<-  validation  region  at  t,  is  defined  as 

k 

I  •  _  i 

.  k  -  j*  :  I*  -  *<tk  tfc.j))  S  (tk  tkl)  [z 


where  d  la  taken  from  the  chl-aquare  table*  and 


klVi'1  i“J 

(4.5) 


■<‘k  Vi’  •  H  “VVi’  ’•"k-i  'k-,’ 


(4.6) 


1 1 


r> .  l  i  ra  I  t  «-d  Memory  Procedures 

Since  the  full-memorv  Probabilistic  Data  Association  Filter  described  in 

N 

Sections  )  and  4  has  growing  memory  and  computational  requirements  (2  hypothe¬ 
ses  can  be  generated  after  N  detections,  each  oi  them  requiring  a  Kalman  filter) 
it  is  not  feasible  for  implementation.  This  leads  to  the  investigation  of 
suboptimal  algorithms. 

One  suboptimal  tracking  scheme,  which  eliminates  the  exponentially  in¬ 
creasing  memory  requirements  of  the  full  memory  filter  can  be  obtained  along 
the  lines  of  the  PDAF  [B1,B2|,  which  has  also  been  called  "zero-st an-back" 
filter  [SI].  In  this  ease  instead  of  all  the  hypothesis-conditioned  estimates 
(i.e.  conditioned  on  a  specific  past  history),  a  single  estimate  is  generated  and 
retained.  Thus,  unlike  the  full  memory  procedure,  in  this  approach,  to  be  called 
"combined  single  hypothesis”  (CSH),  one  does  not  have  to  store  track  history 
hypothesis.  In  this  procedure  when  a  detection  is  obtained  the  probability  of 
the  corresponding  measurement  being  correct  is  computed  and  the  state  estimate 
and  its  covariance  are  updated  using  this  information. 

The  approximate  sufficient  statistic  that  replaces  the  past  data  in  the 
PDA  consists  in  this  case  of  the  following: 


(1)  the  predicted  state  xft^  t^j)  anc*  associated  covariance  matrix; 


(li)  the  arrival  times  8 

k 

t  . 
k 


J»0,1  of  the  latest  detection  of  tvpe  j  prior  to 


In  the  full  memory  approach  these  times  are  obtained  from  the  history  hypo¬ 


theses.  In  the  CSH  procedure  the  arrival  time  of  the  last  t argot -or i r  inat ed  de¬ 


tection  prior  to  time  k  will  Instead  be  determined  from  the  probabilistic  data 


association  as  follows.  The  last  time  the  probabil ist ic  data  association  cal¬ 


culated  that  the  probability  of  a  detection  having  originated  Indeed  from  tlu- 
target  was  greater  than  the  probability  It  was  a  false  alarm  will  be  « I «  «  lured 
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to  be  d,  That  is,  the  most  recent  time  that  a  detection  was  more  likelv  to 

k,  1 

have  been  associated  with  the  target  than  with  background  clutter  is  considered 
to  be  the  last  occurrence  of  a  target-originated  detection: 


d,  ,  ■  t  where 
k,  1  i . 


1  ■  max 

i*k 


(5.1) 


Similarly,  0,  ,  the  last  occurrence  prior  to  t  of  a  false  detection  is  taken 

k.O  k 


as 


0.  ,,  “  i  where  t  -  max  (i  :  6.  n  >  P.  ,  I 

k.o  .0  o  1<k  1.0  i.l 


(5.2) 


When  a  detection  occurs  two  hypothesi s-cond 1 t ioned  estimates  are  calculated. 
The  conditioning  is  an  origin  hypothesis  of  only  the  current  detection  which 
yields 


;<tk  v-kV 


x<Vtk-i) 


(5.3) 


and  x(t,  t,  ,x.  ,)  obtained  from  a  standard  Kalman  filter, 
k  k  k,l 

The  final  state  estimate  is  then  a  weighted  sum  of  these  two  hypothesis- 

ond it ioned  estimates.  The  weighted  factors  are  the  a  posteriori  probabilities 

that  the  origin  hypothesis  is  correct. 

1 


x(t,  *  2-w  *(t.  t..X,  .)  »,  . 

k  k  ~  k  k  k.l  k.J 


(5.4) 


where 


“k.J  '  ’’‘Vj'V'V-  i(,k  lk-l  ’•  •k.o-'k.l* 


3-0.1.  (5.5) 


Thus,  the  probabilistic  data  association  of  the  CSH  algorithm  calculates  the 
probability  of  an  origin  hypothesis  of  only  the  current  detection  being  correct. 
In  the  full  memory  procedure  the  probability  of  a  hypothesis  as  to  the  origin  of 
each  detection  from  initial  to  present  time  must  he  determined.  The  two  pro¬ 
babilities  of  (5.4)  are  calculated  using  Bayes'  rule  as  in  Section  3. 


1  i 


This  completes  the  description  of  the  probabilistic  data  association  for 
the  CSH  procedure.  The  covariance  of  the  state  estimate  (5.4)  must  be  deter¬ 
mined.  Equation  (4.1)  can  be  used.  However,  due  to  the  nature  of  the  CSH 
approach,  the  following  simpler  expression  for  the  covariance  of  the  uncondi¬ 
tional  estimate  is  obtained. 


F<tk  V  ■  \.o  p<tk'Wo>  *  V.  P(tkitk-<k.i> 


*  rk.O  rk,l  uk.l  Uk,l 


(S.b) 


where 


k.l 


x(t 


V*k,l 


)  - 


*(tk!tk-l 


(5.7) 


Beyond  this  simplest  approach  there  are  a  number  of  more  complex  filters  to 
he  discussed  next. 

A  second  subopt lmal  filter,  which  was  Imp  1 ement ed ,  Is  one  In  which  a  fixed 
maximum  number  N  of  hypotheses  are  retained.  For  example,  once  the  number  of 
hypotheses  exceeds  N,  onlv  the  N  most  likelv  ones  are  retained  and  the  others 
are  discarded.  The  probabilities  of  the  retained  hypotheses  are  then  renormal¬ 
ized  such  that  thev  sum  up  to  unity.  Such  a  scheme  is  simple  to  implement  and 
will  give  the  optimum  performance  when  the  retained  probability  mass  is  near 
unity. 

Other  approaches  can  also  be  taken:  one  can  keep  a  predetermined  probability 
mass,  e.g.  99t,  regardless  of  the  number  of  hypotheses  to  be  retained  or  one  can 
combine  hypotheses  that  are  similar  by  some  relevant  standard. 
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6 .  Simulation  Results 

The  problem  which  was  simulated  was  the  tracking  of  a  target  with  a 
four-dimensional  state,  consisting  of  position  and  velocity  in  two  cartesian 
coordinates,  driven  by  white  noise. 

Using  the  notations  of  eqs .  (2.1)  -  (2.4),  the  system  is  characterized 
by  the  following  state  vector 

x  ■  I  t  C  n  n  ]  '  (6.1) 

where  £  and  n  denote  the  cartesian  coordinates  and  the  system  matrices  are 


‘o  1  0  o' 
0  0  0  0 
0  0  0  1 
0  0  0  0 


(6.2) 


B 


0  0 
1  0 
0  0 
0  1 


(6.  3) 


The  process  noise  covariance  matrix  was  taken  as 


Q  - 


(6.4) 


and  the  measurement  matrix 


(6.5) 
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with  measurement  noise  covariance 


9  x  10 


9  x  10 


(6.6) 


The  probability  densities  of  the  inter-arrival  times  of  both  the 
correct  measurements  as  well  as  the  false  alarms  were  chosen  from  the 
gamma  density  family 


f,(t)  XT  (X  t)®'1  e'  V 

(m-1).’ 


t  >  0  (6.7) 


for  the  target  -  originated  detections  and  for  th*  false  detections 


fo(t)  "  (tTT:  (*fVs  °  1  e  F  *  t  >  0 


(6.8) 


where  X_  Is  In  units  (time)  X  is  In  units  ( t lme . volume)  *  and  V  is 
T  r  s 

the  size  of  the  surveillance  region.  The  average  time  between  detections 

originating  from  the  target  Is  m/X  from  the  false  alarms  Is  ®/ ( X  V  )  . 

T  F  s 

The  following  were  the  parameter  choices: 


m  ■  2 


(6.9) 


(6.10) 


and  Xp  assumed  a  number  of  values  to  simulate  various  degrees  of  environmental 
Inter  f  erence . 


The  validation  region  size  was  set  such  that  th*  probability  that  a 


detection  originating  from  the  target  falls  In  It  was  .999. 
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The  threshold  In  (4.5)  was  set  assuming  a  chi-square  densitv,  as 

d  -  13.8  (6.11) 

Three  filters  were  Implemented: 

(i)  a  standard  Kalman  filter  (SKK)  which  processed  any  detection 

lying  within  Its  validation  region  as  if  it  'ere  target  originated; 
(11)  the  CSH  filter 

(111)  a  fixed  memory  filter  In  which  the  N  most  likely  tracks  were  retained. 

The  purpose  of  Implementing  the  standard  filter  is  to  compare  Its 
performance  against  the  two  PDA  filters. 


The  comparison  was  performed  via  a  Monte  Carlo  simulation  for  various 
clutter  densities.  Each  run  lasted  30  units  of  time  and  an  average  was 
obtained  over  100  runs.  Measures  used  to  compare  the  performance  of  the  filters 
Included  the  Cartesian  norm  of  the  position  and  velocity  estimation 
error  at  time  t^  of  the  last  detection. 


H  (x  (t  ) 

p  L 


x(tLjtL)r 


(6.12) 


!  I*vl  I  ■  I  I  Hy[*(tL)  -  *(tL|tL)i|  | 


(6.13) 


where 


H 

P 


10  0  0 
0  0  10 


(6.14) 
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H  .[°100' 

v  [o  0  0  1_ 

The  normalized  position  and  velocity  errors 


(6.15) 


c 

P 


!’"l>  -  1Kp  P  <*!>,.) 


V 


-1 


U.16) 


and  t  with  a  similar  expression  were  also  evaluated.  This  was  done  to  check 

v 

the  consistency  of  the  filter,  i.e.,  whether  its  errors  were  compatible  with  the 
calculated  covariances.  The  expected  values  of  e  and  are  2  because  they 
have,  theoretically,  a  chi  square  distribution  with  two  degrees  of  freedom, 
assuming  no  lost  tracks. 

In  addition,  the  number  a  times  the  filter  "lost"  the  target  was 
determined.  A  target  was  coneidered  to  be  lost  when  more  that  25Z  of  the 
actual  target-originated  detections  fell  outside  a  filter's  validation  region. 

The  parameter  of  our  study  was  the  denalty  of  false  alarms.  The 
parameter  entering  into  (6.8),  X^.  ,  was  taken  as  5,  10,  20  and  30 

for  -  1600.  To  get  an  idea  of  the  severity  of  these  environments  Table  6.1 
shova  the  average  number  of  false  detections  in  the  validation  region  of  the 
standard  filter. 


XF  V. 

5 

10 

20 

30 

Average  Number  of 
falae  detectlona 

1.3 

2.7 

5.2 

7.3 

Expected  number  of 
correct  detections 

15 

15 

15 

15 

Table  6.1  Effective  rate  of  false  va.  true  detectlona 
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The  percentage  of  lost  tracks  Is  presented  In  Table  f>.2.  As  It  can  be 
seen,  for  a  relatively  low  rate  of  false  alarms  the  SKK  has  already  a  signifi¬ 
cant  probability  of  losing  the  track.  The  number  of  lost  tracks  for  both  the 
CSH  and  fixed  memory  filters  is  substantially  lower  than  the  SKK  for  all  clutter 
dens l t ies . 


Filter  v. 

5 

10 

20 

30 

SKF 

.17 

.31 

.64 

.77 

CSH 

.02 

.05 

.09 

.19 

Fixed  eeeory  N  -  8 

.02 

.04 

.07 

.11 

Table  _6._2  Probability  of  Losing  the  Target 


Tables  6.3  and  6.4  present  the  posit  ion  and  velocity  errors,  respectively, 
for  the  three  filters.  For  each  clutter  density,  in  both  position  and  velocity, 
the  SKF  has  the  largest  error  Cartesian  norm  followed  by  the  CSH,  with  the 
fixed  memory  filter  having  the  lowest  error  norm  in  each  case.  In  addition,  the 
SKF  is  very  "optimistic":  the  real  time  covariances  are  much  too  small.  The 
performance  of  the  standard  filter  is  degraded  and  the  covariance  of  the 
estimate  does  not  increase  to  correspond  to  this  situation.  The  normalized  errors 
of  both  PDA  filters  are  larger  than  the  expected  value  of  2  (based  on  no  lost 
tracks)  which  is  due  to  the  times  these  filters  actually  lost  the  target. 
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The  two  PDA  filters  discussed  above  use  both  locations  and  arrival  times 
to  estimate  the  state  of  the  system.  To  determine  the  effect  of  using  arrival 
time  information  on  tracking  performance,  a  "combined  single  hypothesis  with  no 
time"  filter  (CSH-NT)  was  run  incorporating  only  location  information  In  the  PDA. 

The  comparison  of  tracking  ability  for  the  two  CSH  filters  is  displayed  in  Table  6.5. 


1  s 

Filter 

5 

10 

20 

30 

SKF 

.17 

.  31 

.64 

.77 

CSH-NT 

o 

_ 

.18 

.  33 

.52 

CSH 

.02 

.05 

.09 

.  19 

Table  6.5  Probability  of  Losing  the  Target 


These  results 
realized  from 


indicate  that  significant  Improvement 
use  of  arrival  time  Information. 


in  tracking  ability  can  be 
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7 .  Conclusion 

An  algorithm  has  been  developed  to  track  a  target  where  the  origin  of  any 
detection  la  uncertain.  In  addition,  the  detections  do  not  occur  at  fixed  Intervals, 
but  Instead  the  time  Intervals  between  detect  Iona  are  random.  Since  the  optimal 
procedure  requlrea  exponentially  Increasing  memory  and  computation, auboptlmal 
versions  are  presented.  A  simulation  was  performed  to  compare  the  performance 
of  these  algorithms  versus  the  standard  filter. 

The  least  expensive  filter  that  uses  a  PDA  procedure  -  the  combined  single 
hypothesis  (CSH)  algor lthm  -  provides  a  substantial  performance  improvement  over 
the  standard  filter.  The  fixed  memory  filter,  which  carries  the  N-8  most  likely 
hypotheses,  yielded  a  performance  which  is  probably  verv  close  to  optimum.  It 
provided  further  Improvement  over  the  CSH,  but  a'  a  substantial  cost  -  it  Is 
8  times  more  expensive  than  the  CSH.  The  trade-off  between  performance  and  cost 
Is  a  question  which  must  be  Judged  within  the  context  of  an  individual  problem. 

The  deletion  of  the  time  of  arrival  Information  from  the  CSH  filter  resulted  In 
significant  performance  deterioration.  This  Indicates  that  the  time  of  arrival 
Information  can  be  used  successfully  to  obtain  significant  tracking  performance 
Improvement  In  this  class  of  problems. 
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1.  Introduction 

The  problem  of  associating  probabilistically  measurements  to  a  number  of 
known  targets  (established  tracks)  was  first  considered  in  ( B 1  ] .  The  resulting 
algorithm,  which  was  an  extension  of  the  results  of  [B2],  used  a  diffuse  prior 
for  the  false  measurements.  While  this  was  shown  to  be  equivalent  in  the 
single  target  case  to  a  Poisson  model  (see  [ B 3  ] )  it  was  not  adequate  for  the 
multiple  target  case.  More  recently,  an  algorithm  was  presented  In  ( R 1 )  that 
yielded  probabilities  of  measurement  histories  being  a  track.  While  this 
latter  algorithm  Is  more  general  than  all  the  previous  ones,  it  has  a  number 
of  limitations.  Its  feasibility  of  application  for  complex  state  models 
where  the  state  is  not  completely  observable  from  one  sensor  is  an  open 
question.  The  model  used  for  "new  targets",  which  consisted  of  an  arrival 
process  with  strict  lv  positive  rate  can  lead  to  saturation  of  the  surveillance 
region. 

The  algorithm  to  be  presented  in  Section  2  deals  with  the  association  of 
an  arbitrary  number  of  measurements  to  an  arbitrary  number  of  established 
tracks  (targets). 


T 


For  the  sake  of  Illustration  consider  first  the  case  with  n-2  known 
targets. 

Let  the  validation  regions  ("gates")  of  the  two  targets  be 

Vj  -  U  :  ||  z  -  ||  <  y)  J-1,2  (1) 

These  regions  are  depicted  In  Fig.  1  with  a  total  number  of  validated  measure¬ 
ments  10*4  Including  one  which  Is  common. 


/ 


Fig.  1  Two  validation  regions  with  nonempty  Intersection 
Next,  consider  the  general  case. 

The  validation  matrix  for  m  measurements  and  n  targets  is 

$)  -  [a^  ]  1-1 . m  J-0,1 . n  (2) 

with  binary  elements  Indicating  whether  measurement  1  has  been  validated  for 
target  J.  Index  J-0  stands  for  "clutter".  The  validation  matrix  corresponding 


to  the  above  figure  Is  (for  n-2) 
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(3) 


The  probabilities  of  Joint  events 


-n 

1-1 


(«> 


where  *  means  that  measurement  1  originated  from  target  J  ,  will  be  computed. 
J1  1 

To  obtain  the  feasible  events  we  shall  use  an  Indicator  variable  equivalent  of 


Xtj .  namely 


|  1  If  Xy  occurs 
I  0  otherwise 


(5) 


Then  to  a  joint  event  ^  there  corresponds  a  matrix 

ft(*)  -  1  (6) 

The  conditions  for  an  event  to  be  feasible  are 

(1)  only  one  unit  from  each  row  of  Q  can  be  picked  up  for  0 


l  (l)  ’  1  1*1.  •  •  •  .®  (7) 

J-0  J 


(11)  at  most  one  measurement  can  belong  to  each  target 

,  m 

6 a  <A>  *  I  w..  (l)  <  1  J-l , . . . ,n  (8) 

J  1-1  lJ 


The  binary  variable  5^  Indicates  whether  target  J  was  detected  and  Its 
measurement  validated  In  event 
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Measurement  1  Is  associated  with  a  target  In  event  ^  if 

v*>  ■  l  “11  <*>  -  1  <9> 

1  j-i 

An  algorithm  that  obtains  all  the  feasible  matrices  1)  can  be  set  up  as 


follows 


1.  Scan  through  the  (say,  N)  non-zero  elements  of  Q  to  set  up  all  the 
possible  binary  numbers  (2‘  )  from  them 

2.  Eliminate  all  those  that  do  not  satisfy  (7)  and  (8). 

Using  ft  one  can  obtain  the  set 


Xlj  ‘  <  A  :  “ij  (*)  "  1  } 


needed  to  evaluate  the  marginal  probability 


8tJ  -  p  (xtJ!z 


ft)  -  X*  aI zk.  n> 

*LXn 


from  the  probabilities  of  the  Joint  events 

The  probabilities  of  the  events  ^  we  obtained  from 

P<Alzk,u)-  j  p(z(k)|A,zk"1,n)  P{nU.zk_1}  p (jL|zk_1l  <«) 

k  k-1 

where  Z(k)  the  set  of  current  observations,  Z  “  Z(k)UZ  is  the  set  of 
observations  through  k. 

The  first  term  on  the  r.h.s.  of  (12)  is  obtained  next 


plz(k)|A.zk‘1,nj  -  |-j  p l*t (k) I xtJ .Z11*1  ,n) 


where 


k-1  I  C1r  fil*l(k)l 

Pl*1(k)|xlj.z  1,H]  -  <  J  1  J 

(  v‘l 

\  r. 
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In  the  above  f  (•)  *8  t*'e  density  of  the  predicted  measurement  for  target  J  (at 
.he  current  time  k) ;  c  *  f.  is  this  density  restricted  to  the  region  F  in 

Jri  J  ri 

which  z ^ (1c)  lies.  For  example  in  Fig.  1  the  three  regions  V  are  illustrated. 
There  is  one  such  region  for  each  different  row  of  SI  and  their  total  number 
is  ra^.  For  the  validation  matrix  (3)  the  indices  r^  corresponding  to  each 
measurement  in  Fig.  1  are 


3 


The  false  measurements  pdf  in  (14)  is  the  Inverse  of  the  volume  of  the 
corresponding  region 


>1  U'  l 
I  M 


(13) 


The  constralts  c  are  given  by  the  expression 
lri 


Cjr.  ‘  /  V' 


jl-l(k)]  dt1(k) 


(16) 


ri 


which  does  not  have  to  be  calculated  because  it  will  cancel  out  as  shown 
later. 

The  second  term  on  the  r.h.s.  of  (12)  is 

P{S%,Zk-1>  -  P  n  U.(k)cP  |Xl1.Zk_1> 

Ll-1  1  1  J 


•(I  I 


p(«1(k)«Vri|xljf* 


k- 


>))p[  o 


(«1(k)cpf. 


l*ioZ 


k-1, 


(17) 


where  the  product  is  over  the  measurements  associated  with  a  target  (since  the 
corresponding  events  are  Independent);  the  second  term  contains  the  measure¬ 
ments  not  associated  with  any  target. 
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Note  that  for  a  detection  associated  with  a  target  the  probability  that 
the  corresponding  measurement  falls  in  V  ,  conditioned  on  Xji  Is 

ri  r  J 

Pl*t(k)e^r  Uij.Z11"1}  -  /  fJl*1(k.)]d«1(k>  /  fjlz1(k)]dzi(k) 


if  ^(ji)  "  1 


The  denominator  in  the  above  follows  from  the  fact  that  if  is  considered 

then  z^  was  validated  for  target  J  and  thus  its  pdf  is  f  divided  bv  the 
probability  that  it  fell  in  the  gate 


p  4  /  fj(«> 

J  V. 


Then  using  (16),  Eq.(18)  becomes 

P<«1<k>«*rtlxiJ.2k'l>  ‘  <=lrj  ^  <20> 

The  last  term  in  (17)  can  be  rewritten  as  the  product  of  probabilities 

that  the  number  of  false  measurements  if  in  the  nonover lapping  regions  V  , 

ri 

originating  from  a  Poisson  distribution  with  rate  X,  are  as  specified  by 
the  event  £  under  consideration 
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The  last  term  from  (12)  Is  the  prior  (w.r.t.  the  current  time  k) 
probability  of  event  Note  that  (17)  Is  the  probability  that 

(I)  the  correct  measurements  Implied  by  ^  (l.e.  given  that  they  are 
detected  and  validated)  fell  into  the  corresponding  nonoverlapping 
regions  T 

ri 

(II)  the  incorrect  measurements  fell  Into  the  various  regions  V  as 

Implied  by  In  volume  ,  l.e.,  their  number  and  apportionment 

among  the  regions. 

In  view  of  this  the  last  term  will  reflect 

(1)  the  probability  of  detection  and  validation  of  the  correct 

measurements 

(2)  the  probability  that  the  measurements  of  the  targets  whose  6^-0 
(no  measurement  associated  with  target  J)  were  either  not  detected 
or  fell  outside  the  validation  region. 

It  will  be  assumed  that  all  Joint  events  that  have  measurements  from  the 
same  set  of  targets  are  a  priori  equiprobable. 

Denote  the  binary  vector  that  Indicates  which  targets  are  assumed 
detected  in  an  event  ^  as 

6  -  6  Cl)  (23) 

with  components  defined  by  (8).  There  are  2n  such  different  vectors  where 
n  is  the  number  of  targets.  All  Joint  events  £  that  have  the  same  £  arc 
thus  considered  a  priori  equiprobable. 

Let 


qfi  "  q(£,n) 


(24) 


be  the  number  of  events  ^  with  the  same  &_  as  It  follows  from  the  matrix  ft  (see 
Appendix  A  for  the  expression  of  this  number  and  consistency  check  with  [ B2 ] ) . 
With  this  the  prior  probability  of  an  event  ^  is 

PU|zk-‘.  -i-  I  I  p  p  I  I  (i-,  r  )  <»> 


Note  that 


E  Pli|zk'1' '  £  r  I  V  pc.  r  IK  pc.) 

A  i  J  J  j:6j(i)-o  3  3 

•  n  (\  \  *  -o,  s )  • p 


which  indicates  that  the  normalization  constants  in  (25)  are  "natural". 
Combining  (14),  (22)  and  (25)  into  (12)  yields 


P{A|Z* 


•al  ‘  \  I  I  V,  V*t' 


1:t1(2.)“1 


i:T1(ji.)-o 


v  -IV, 


1  K/c!  [I'  ‘  <»V 


*  n,!)-1 


t,  1  1  P»,  S  1  1 


(‘•'-J  S) 


y\  <*. 


•  I  I  p»,  I  I  K  *<■') 
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where 


v 

♦  \ 

t-1 


(A) 


and  c  ,  V  ,  P  cancel  lmaedlately  and  the  tenn 

iri  GJ 


n  «'1Vt 


.x*. 


-w-  -A  t-l  -XV 

e  -  e 


t-1 


(28) 


(29) 


appears  regardless  of  which  ^  la  under  consideration  and  thus  cancels  between 
numerator  and  denominator  In  (27). 

The  probability  of  a  Joint  event  ^  given  in  (27)  was  derived  assuming 
that  there  are  validation  regions  corresponding  to  each  target.  While  the 
volumes  of  the  validation  regions  cancelled  in  the  final  equation  some 
factorials  corresponding  to  the  number  of  false  detections  subsumed  in  ^ 
were  left  in  the  final  expression.  Also  the  normalization  constant  defined 
in  (24)  appears  in  (27). 

The  derivation  of  the  conditional  probability  of  a  Joint  event  is  done 
next  assuming  each  validation  region  to  coincide  with  the  entire  surveillance 
region.  In  this  case  the  validation  matrix  (2)  is  made  up  of  units  everywhere. 
(However,  for  practical  purposes  one  can  still  use  validation  regions  to 
avoid  considering  events  that  will  have  negligible  probability.)  In  this 
case  there  is  no  ft  in  (12),  which  becomes 

p{iiZk>  -  i  p[Z(k)|*,Zk-1]  P(i|Zk_1)  (30) 

The  first  tens  above  is 

p(Z(k)|*,Zk_l]  -  pj  pl«t(k)  | Xtj •  Zk— 1 )  (31) 

1-1 


where 
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ki  Kim*)) 
pl^Wlx^.z  i  -  3 

(  V-1 


(32) 


with  V  being  the  volume  of  the  entire  surveillance  region. 
The  "prior"  probability  of  an  event  ^  is 
.k-] 


H*|z 


'  $3  1  I  p0j  I  I  ('-%) 

J:6  (*J-i  3  J:«j(A>-o  3 


_XFV  4  -1 

•e  F  (XpV)*  (♦!)  1 


(33) 


where 


♦  -  ♦(*.)  (3^>) 

is  the  total  number  of  false  measurements  in  event  ^  and  a(j^)  is  the  number 
of  events  ^  with  detections  from  the  same  set  of  targets.  This  is  equal  to 
the  number  of  permutations  of  m  (total  number  of  detections)  takes  as  m-i 
(number  of  correct  detections) 

(35) 


a(*>  “  C#  ‘  “  fr 


Inserting  ( 31 )  —  ( 3 5 )  into  (30)  yields 

p‘il2k  1  -  £  I  I  ft‘*£ 
1:^(1)-! 


(k)  ]  V 


■  a  n  v  i  i  (i-\ ) 

J:«j(i)-1  3  JsSj^-o  J 


■XFV  A  -1 

.•  *  (xFv)*  (♦:) 


(36) 


which  becomes  after  cancellations 


n 

1:^(2.)- 
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3.  An  Example 


Consider  two  targets  with  a  one-dlmenslonal  measurement  with  predicted 
locations  for  targets  1  and  2  as 

*l"°.  *2"3 

and  variance  of  the  prediction  S“l.  The  following  three  measurements  were 
to  be  associated  with  the  two  targets 

r  ■  -1.5,  i,  ■  1.5,  i  -4.5 
a  b  o 

This  Is  depicted  below.  Only  z^  was  validated  for  both  targets. 


Table  1  presents  the  probabilities  of  all  the  Joint  events,  assuming 
Pp  •  0.8,  for  various  densities  X  for  the  false  alarms.  Table  2  shows  the 
corresponding  marginal  probabilities. 

Note  that  in  the  absence  of  target  2,  the  events  xa  ^  and  xb  j  would 
have  been  equlprobab le .  Since  i.  can  also  belong  to  target  2,  the  Joint 

D 

PDA  yields  a  drastically  different  result.  This  Illustrates  the  need  to 
account  for  all  the  possible  origins  of  each  measurement. 


X-.05  I  X-.02 


.0284 


70004 


7945  I  .3168 


Table  1.  Joint  probabilities  for  all  the  feasible  associations 
for  various  false  alarm  densities. 


Cone  lug  Ion 


Raid  [Rl]  derived  his  expressions  for  the  joint  probabilities  assuming 

(I)  The  target  detection  events  are  Identically  distributed  (same 
Pp)  and  independent 

(II)  A  number  of  new  targets  are  detected  according  to  a  Poisson 
distribution  (this  Implies  a  net  Inflow  into  the  surveillance 
region,  l.e.,  monotomically  Increasing  number  of  targets) 

(III)  Initialization  of  new  targets  from  one  measurement  is  possible 
The  JPDA  approach  described  above  Is 

1.  More  general  than  Reid's  because  the  target  detection  events  are 
not  Identically  distributed 

2.  Less  general  because  It  assumes  no  new  targets 

3.  More  realistic  because  both  assumptions  (11)  and  (111)  above  do  not 
seem  to  fit  many  practical  situations,  in  particular  the  passive 


sonar  tracking. 


A-l 


Appendix  A 


In  the  particular  case  of  n-1  target  the  computation  of  the  priors  Is 
done  as  follows.  For  x^.  l“l,...,m  ("the  1-the  validated  measurement  is 
correct")  one  has  6-1  ("target  detected")  and  from  (24) 


qj  -  m 


(A.  1) 


From  (25)  It  follows  chat 


PlXl|2k-1)  -iPDPc 


(A. 2) 


For  l-o  one  has  6-o  and  q  "1,  thus 

o 

p(x0|zk_1}  -  i-pd  pg 


(A.  3) 


Computation  of  the  posterior  probabilities  Is  done  next.  Note  that  In  (14) 
f  Is  the  nontruncated  density  of  the  predicted  measurement.  In  the  survey 
paper  [B3 ]  f  was  used  as  the  truncated  density,  l.e.,  it  included  the  factor 

P,1.  From  (27)  one  has  for  it*o  (In  which  case  ♦-m-1,  q  -m) 

G  o 


k  \  *  1 

6l-P{Xl|7.  .  n^-m)  -  f(*1)  -^yr  PD 


(A. 4) 


for  l*o 


0o-P(xo|z  .  mk-«)  ■ 


Am  1 

rr  ~T  (1-Pn  Pr) 

C  tB,  U  (.» 


(A. 5) 


The  normalization  constant  c'  Is 


»-l  _A_  i® 

’  -  V  pd£  f(V  (1-pn  PG> 


(A.  6) 


Then,  for  1-1,... ,i 


sr  pp  t(«i>  f(ti> 

1  "  Is  X**1  ^  X=^ 

h  (1-W  +  ^  ^^i5  bo  +  V 


(A.  7) 


A-2 


where 


1-P  P 
D  C 


(A.  8) 


v 

t 


f(zj 


(A. 9) 


which  matches  the  earlier  results  of  [B2], 

Consider  n-2  targets  with  the  foll>iwing  number  of  validated  measurements 
-  for  target  i  only  i-1,2 

m^j  -  for  targets  1  and  J  (in  the  Intersection  of  their  validation 
regions) 

Then  the  number  of  events  for  different  _6  ■  (6^  6^]  is 

q(5l-l,  «2-l)  -  mu  »22  ♦  -n  ^2  ♦  »12 

q(*1"o.  42“1)  “  “12  +  “22 

q(6l“l-  Vo)  •  “11  +  “12 

q(4j-o,  «2-o)  -  1 
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